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Abstract — This paper derives recursion equations for a robust 
smoothing problem for a class of nonlinear systems with uncer- 
tainties in modeling and exogenous noise sources. The systems 
considered operate in discrete-time and the uncertainties are 
modeled in terms of a sum quadratic constraint. The robust 
smoothing problem is solved in terms of a forward-time and a 
reverse-time filter. Both these filters are formulated in terms of 
set-valued state estimators and are recast into subsidiary optimal 
control problems. These optimal control problems are described 
in terms of discrete-time Hamilton-Jacobi-Bellman equations, 
whose approximate solutions lead to recursive Riccati difference 
equations, filter state equations, and level shift scalar equations 
for the forward-time and the reverse-time filters. 

I. Introduction 

Filtering and smoothing are extensively used for estimation 
of data from noisy signals. Although filtering by itself is 
sufficient in most cases, smoothing becomes necessary where 
outliers need to be estimated from noisy data. These topics 
have been of interest in the estimation and control community 
since the 1960s; see e.g., 0, 0, 0, 0. In spite of the 
advancements in the theory for estimation and smoothing, 
the success of most of the algorithms are based on precise 
knowledge of the system dynamics. In addition, such methods 
are known to diverge in the presence of unmodeled system 
dynamics and parametric uncertainties 0. In order to take 
into account uncertainties in modeling, robust versions of 
these estimators were developed in the stochastic as well as 
deterministic settings; see e.g., 0, Q, 0, 0, HO), E], 
Ifl2l, Ifl3l . HI, |Q3J, H6l, ED. Also, the issue of robust 
smoothing for systems with model uncertainties have been 
addressed in the stochastic and Hoo frameworks in [18|, [19|, 

lEol, ED, ED, ED. 

As seen from research in the past, most of the robust 
filtering, prediction, and smoothing techniques are not ap- 
plicable to systems where the statistical nature of noise is 
unclear and hence an appropriate statistical model for the 
system is not available. Such a situation may arise in systems 
where the uncertainties in modeling dominate over uncer- 
tainties due to exogenous noise sources. In such cases, a 
deterministic approach to filtering and smoothing might be 
more appropriate. Although there are various approaches to 
deterministic filtering, [24] provides a set-membership ap- 



proach to Kalman filtering. This set-membership approach 
to deterministic filtering involves computing the set of all 
possible states compatible with given output measurements 
for an uncertain system with norm bound uncertainties. This 
method was extended by Savkin and Petersen in 0, |8| to 
linear uncertain systems where the uncertainties were bound 
by integral and sum quadratic constraints. Furthermore, the 
method in 0, was extended to include nonlinear uncertain 
systems in 0, El, 03J, H6), ED. However, all these 
methods deal with the problem of deterministic filtering for 
uncertain systems and do not deal with issues related to robust 
smoothing. The set-membership state estimation approach in 
0, was used to derive a robust smoothing algorithm for 
uncertain linear systems in [25 1. This method was extended 
to uncertain nonlinear systems in the continuous-time case in 

E3- 

In this paper, we present a robust smoothing algorithm for 
uncertain nonlinear systems in a discrete-time setting, with 
uncertainties bound by a sum quadratic constraint (SQC). Such 
a smoothing algorithm is applicable for systems in fields such 
as aerospace engineering and signal processing where discrete- 
time sensors are used. In this paper, we derive a fixed-point 
smoothing algorithm that solves a bridge filtering problem. 
This bridging approach to smoothing comprises of solving 
two underlying robust filtering problems. The corresponding 
solution leads to a forward-time filter over the discrete-time 
interval [0, k] and a reverse-time filter over the discrete-time 
interval [k,i\. Here, < k < t. 

Both the forward-time and reverse-time filtering problems 
are formulated as set-valued state estimation problems which 
are then recast into optimal control problems. The optimal 
control problems are described in terms of Hamilton-Jacobi- 
Bellman (HJB) equations evaluated over different time in- 
tervals. Various nonlinear systems in the HJB equations are 
linearized and approximate quadratic solutions are computed 
for the dynamic programming equations. This leads to a 
Riccati equation, a filter state equation, and a level shift scalar 
equation for both the forward-time and the reverse-time filters. 
Although it is straightforward to formulate the set-valued state 
estimation problem for the reverse-time filter, this is not trivial 
in the case of the forward-time filter. The set-valued state 



estimation problem for the forward-time filter will be solved 
by obtaining an approximate solution to the optimal control 
problem in terms of a forward-time dynamic programming 
equation for a reverse-time discrete-time uncertain nonlinear 
system with uncertainties described in terms of an SQC. For 
a detailed explanation, see ll27ll . 

The rest of the paper is organized as follows: Section 
ITT1 introduces the formulation of the discrete-time uncertain 
nonlinear system starting from an equivalent continuous-time 
system. The set-valued state estimation approach to determin- 
istic filtering is presented in Section [Ell] The formulation of 
the robust smoothing problem in terms of a forward-time and a 
reverse-time filter are presented in Section |IV] The recursion 
equations for the discrete-time robust nonlinear smoother in 
terms of the Riccati equation, the filter state equation, and the 
level shift scalar equation are derived in Section |V] The paper 
concludes with a summary and future research directions in 
Section EB 

II. Problem Formulation 

In this section, we introduce a general model for a 
continuous-time nonlinear uncertain system where the un- 
certainties are modeled in terms of an integral quadratic 
constraint (IQC). Furthermore, we provide a discussion on 
obtaining an equivalent discrete-time uncertain system where 
the uncertainties are modeled in terms of an SQC. 

Consider the following forward-time continuous-time un- 
certain nonlinear system 

x(s) = a c (s, x(s), u(s)) + D c w(s), 
z(s) = 9c(s,x(s),u(s)), 

y(s) = c c (s,x(s),u(s)) +v(s). (1) 

Here, x(s) G M" is the state vector, z(s) G R 9 is the uncer- 
tainty output, y(s) G M 1 is the measured output, w(s) G M. p 
and v(s) G M 1 are the process and measurement uncertainty 
inputs. Also, a c (-), g c ('), andc c (-) are given sufficiently smooth 
nonlinear functions of the state variables and D c is a matrix 
function of time. Such a system is depicted in Fig. [JJ 
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Fig. 1. Uncertain system. 

The uncertainty inputs w(s) and v(s) are assumed to be 
square integrable, satisfying an IQC at any time t > as in 



0, ED, 

S c (t,x(t),w([0,t}) := \\x(0) - x \\% 



(\\w(s)\ 



\Q {s) + \HsW R{s) -\z(s)\')ds<d. (2) 



Here, xq is the initial state value, xq is the nominal initial 
state, and u>[o,t] := w(s)|o< s <t represents the restriction of 
the noise term w(-) over the time interval [0,t]. All vectors 
in (01 and in the rest of the paper are organized as columns, 
unless specified otherwise. Also, w(s) and v(s) are admissible 
uncertainties described by 



v(s) 



(3) 



where ?;(•) is a nonlinear time varying uncertain function and 
x(-) represents the state vector x(s) for all instants of time 
s G [0,t]. The fact that the uncertainty inputs w(s) and v(s) 
depend on x(-) in (0 indicates that they are allowed to depend 
dynamically on the system state. Also, N = N T y is a given 
matrix and Q(s) = Q{s) T >- 0,R(s) = R{s) T ^ are given 
matrix functions of time. 

In order to derive recursion equations for the discrete- 
time robust smoother, we need to convert the continuous- 
time nonlinear system in (Q~|i and the corresponding IQC 
(fJJ into corresponding discrete forms. As mentioned in the 
introduction, the recursion equations for the smoother will be 
obtained by solving for a forward-time filter and a reverse- 
time filter. The derivation of the forward-time filter in the set- 
valued state estimation framework, is most straightforward if 
the continuous-time uncertain system in ((TJ is discretized in 
reverse-time rather than in forward-time. The application of 
higher-order Euler or Runge-Kutta methods to ([T) gives us 
the following reverse-time discrete-time system 



x t = a t (x t +i) - D t w t , 

zt = g t (x t+1 ), 

y t = ct(x t+ i), +vt+i, 



(4) 



where a^,gt, Ct are discrete-time nonlinear functions and D t 
is a given matrix function of time, all expressed at time step 
t. 

Indeed, the solution to the reverse-time filter in the set- 
valued state estimation framework is obtained by discretizing 
the continuous-time uncertain system (HJ forward in time. 
Applying higher-order Euler or Runge-Kutta methods to ([TJ, 
we arrive at the following forward-time discrete-time uncertain 
nonlinear system 

x t +i = a t (x t ) + D t w t , 
Zt+i = K t (x t+ i), 

Vt+i = +v t +i, (5) 

where a t , K t , ft are given nonlinear functions of time and D t 
is a given time-varying matrix at time step t. The nonlinear 
functions in © and (0 include the effect of known control 
signal u t at time step t. However, their explicit dependence on 



the control signal u t will be omitted as the paper deals with 
a smoothing problem. 

In the case of the forward-time filter and the reverse-time 
filter, the uncertainties associated with the systems © and (O 
are modeled in terms of an SQC. This SQC is derived by 
discretizing the IQC in (O using Reimann sum approximation 
to give 

St(xt,W[ ,t]) 



Fo 



Xo\\ 2 N + J2 



2 



< d. 



(6) 

Here, N — N T >- is a given matrix and Q s = y 
0, R s — Rj y are given matrix functions of time. 

Since the recursion equations for the robust smoother are 
solved in terms of a forward-time filter and a reverse-time filter 
in the set-valued state estimation framework, the following 
section introduces the concept of set-valued state estimation. 

III. Set- valued State Estimation 

In the continuous-time case, consider an output yo(-) for any 
s G [0, t] for the system in ([T). The set-valued state estimation 
problem comprises of finding the set 

Zs[x ,yo(-)\ s ,d] (7) 

of all possible states x(s) at time s for the system ([T|i with 
initial conditions and uncertainties defined by the IQC ©, 
consistent with the output yo{-)- Given an output sequence 
yo(-), it follows from the definition of Z s that 

x s e Z s [x ,y (-)\ s ,d] (8) 

if and only if there exists an uncertain input sequence 
w(-) such that, J[x s ,w(-)} < d, where the cost functional 
J[x S) w(-)] is defined in terms of the IQC © as 



J[x s ,w(-)] = inf S c (t,x(t),w(-))] 

w(-) 



(9) 



with yo(t) = c c (x(t)) + v(t), z(t) — g c (x(t)). Now we present 
the formulation of the discrete-time robust smoothing problem 
in the set-valued state estimation framework. 

IV. Robust Smoothing and the Optimal Control 
Problem 

In this section, we consider a discrete-time robust smooth- 
ing problem formulated in the set-valued state estimation 
framework. Here, the robust smoothing problem is considered 
at a fixed point in time k = (t — t) and comprises of 
computing all the states x k £ X k [xo, y[o,k], d] at time k. 
Indeed, the corresponding optimal control problem in this case, 
will depend on the time interval being considered. 

From the definition for the set-valued state estimator in (0- 
©, at time k we have in the discrete-time case 



X k [x ,y [0tt] ,d} — \x k £ K" : inf (s x ( 
L w io,t] V 



Xk, W[Q,fc]J 



+ S 2 (x k ,w [ktt] y\ < d|, 
(10) 



where the infimum is taken over all xr.) and i«[o,t] with the 
boundary condition xt — x k . The value function Si(-) is 
defined over the time interval [0, k] as 



Sx(xk,W[ 0tk ]) 



\x -x \\ 2 N +Y^ 



(11) 



The value function S^O) is defined over the time interval [k, t) 

as 



S 2 (x t7 w [kA ):^Y / (II^HL + INIL - l^l 2 ) • (12) 

s — k 



Note the absence of the contribution of the error in initial 
state in (V2\ . This is because the initial state is known a priori 
over the time interval [k, t]. Indeed, for each x k , Si(x k , vJ[o,t]) 
depends on the sequence W[o,t] over the time interval [0, k], 
whereas, S2(x k , ifrfc.ti) depends on the sequence Wo,*] over 
the time interval [k,t]. Hence, the left hand side in the 
optimization problem in (fT0T > can be split into the sum of 
two optimization problems based on the time interval being 
considered. The first optimization problem depends on the in- 
fimum computed over Wf.-s £ £ 2 [0, k], whereas the second one 
depends on the infimum computed over ui(.) 6 C 2 [k, t]. Hence, 
the definition for X k [xo, y\o,t], d] i n <n~0T > can be restated as 



x k[xo,V[o,t],d] = e 



inf 

u) ( .)G£ 2 [fe,t 



\Si(x k ,W(.^ 
(^S 2 (xk,W(.))j < d.\ 



inf 

Mi ( .)e£ 2 [o,fc] 



(13) 



The optimization problem 



inf Si(xk,wr.))+ inf S 2 (x k ,wt: ) ) (14) 

for the systems and (O satisfying the quadratic constraints 
defined in (flTT l and (fT2l > respectively, defines a nonlinear 
optimal control problem with a sign indefinite quadratic cost 
function. The robust smoothing problem will be solved by 
obtaining an approximate solution to this optimal control 
problem. 

V. Robust Nonlinear Smoother 

In order to obtain recursion equations for the robust non- 
linear smoothing problem, the two optimization problems in 
(fl4l i will be individually solved using methods of dynamic 
programming. This is accomplished by finding approximate 
solutions to the underlying HJB equations. Note that the 
solution to the robust nonlinear smoothing problem consists 
of two filters: one operating in forward-time over the interval 
[0, k] and the other operating in reverse-time over the interval 
[k,i\. The solution to the forward-time filter will comprise 
of approximately solving the HJB equation forward in time 
for the reverse-time discrete-time uncertain nonlinear system 
in (HJi and the uncertainties described by the SQC ( ITTV The 
solution to the reverse-time filter will involve solving for the 



HJB equation reverse in time with the forward-time discrete- 
time uncertain nonlinear system (0 and the uncertainties 
described by the SQC (fT2l . The formulations and solutions 
to the optimal control problems for the forward-time and 
reverse-time filters will be discussed in Sections IV-Al and IV-BI 
respectively. 

A. Forward-Time Filter 

In order to derive the Riccati equation, the filter equation, 
and the level shift scalar equation for the forward-time filter, 
we consider the following optimization problem involving 
Si(-) from (TT~4T > defined over the time interval [0, k], 

H k (x) := inf Sx(x k , Wt.)). (15) 

Here, H k {x) defines a reverse-time optimal control problem 
with Si(-) defined in fTTY In the continuous-time case, this 
optimization problem can be solved using dynamic program- 
ming methods in terms of the forward-time system dynam- 
ics. However, this is not straightforward in the discrete-time 
case, since it would involve inverting the nonlinear function 
describing the discrete-time nonlinear system, which may not 
be possible. In order to solve for the set-valued state estimation 
problem in this case, we consider an approximate solution for 
the optimal control problem (fTBl in terms of a forward-time 
dynamic programming equation for the reverse-time discrete- 
time uncertain nonlinear system (H) with an SQC uncertainty 
description fTTY This technique was first proposed in [14], 

urn 

The corresponding HJB equation for the optimization prob- 
lem in (fTBI l for the reverse-time discrete-time system (2) and 
the SQC fTTJi is given by 

V k+1 (x, y :k) = inf (V k {a k (x) - D k w k ) 

Wr.) G£ 2 [0,k] 

+ IKIIq s + IML-I^I 2 ) 

(16) 

with initial condition 

Vo(x) = \\x-x f N . (17) 

As a first step towards obtaining an approximate solution to 
the HJB fToY we approximate the value function V k ( ) as 

V k {x) « \\x-x k f nk +4>k- (18) 

Applying the approximate solution (fT8l to the HJB fToT ) we 
get 

\\x - x k \\ 2 n +(j) k+ i = inf (\\a k (x) - DhWhlfa 
w() e£ 2 [o,fe] 

+ IKIlQ fc + lkllk-l^l 2 + <^)- (W) 

Solving for the infimum problem on the right hand side of 
(fT~9T > using the method of completing the squares, we obtain 
the following optimum value for w k , denoted as w* 

w* = (Dln k D k + Q fe ) _1 ^n fc (a fe (x) - x k ). (20) 



Substituting w k = w* from fZUl i into (fT~9b and rearranging 
terms we get 

\\ x - x k\\n k+1 +(f>k+i = ||a fc (a;)||| - 2xlY,a k {x) 

+ PfcHl + \\Vk - c k (x)\\ 2 Rk - \g{x)\ 2 + 4>k, 

(21) 

where 



E = n fc - U k D k (DlU k D k + Q k )*DlU k . (22) 



Here, (•)# denotes the Moore-Penrose pseudo-inverse. If the 
matrix (D^H k D k + Q k ) is positive-definite, the pseudo- 
inverse in d22l can be replaced by the normal matrix inverse. 
This is feasible and holds in most cases where a suitable 
solution exists for n^; see e.g., |8]. 

In order to obtain an approximate solution to the right hand 
side of (f2TJ, we consider the following first order linearization 
for various nonlinear functions about the point x k . 

a k (x) w a k + A k (x - x k ), 
g k {x) w g k + G k (x - x k ), 
c k (x) w c fc + C k (x - x k ), (23) 

where a k = a k (x k ),g k = g k (x k ), and c k = c k (x k ). Also, 
A k , G k , and C k are the Jacobian matrices of the maps a k , g k , 
and Cfc respectively. 

Substituting for various linear approximations from (l23l into 
(|2TT > we obtain 

\\ x ~ ^lln fc+1 + <t>k+l = \\ x \\ 2 AlT,A k +ClT,C k -GlG k 

- 2x T ( K {Al^A k + AlV + Cl^C k - G k G k ) x k 

-A T k m k + ClR{y-c k ) + G T k G k ) 

+ \\ X k\\ 2 AlY,A k +T,-2T,A k +ClR k C k -GlG k 

- 2xl(Alm k + Ea fc - ClR k {y k - c k ) - G^g k ) 

+ a^Ea fc + \\y k \\ 2 Rk - 2ylR k c k + ||4||| fc - glh + <f>k- 

(24) 



Comparing terms on the left hand side and right hand side 
of (l24i i. we obtain the following recursion equations for the 



forward-time filter corresponding to the time interval [0, k\: 

Riccati equation over the time interval [0, k] : 

Ilfe+i = AjY*A k + Cj.HCk — G k G k , 

n = N. (25) 

Filter state equation over the time interval [0,k] : 

Xk+i = x k + H, 

xo = xo, (26) 
where 

5 = (AlZ(x k - a fc ) + ClR{y - c fc ) + G T k G k ). 
Level shift scalar equation over the time interval [0, k] : 

4>k+l = <f>k + \\ X k\\\T^ Ak+s _ 2SAk+c T RkCk ^ G T Gk 

- 2xl[Alm k + Sa fc - ClR k {y k - c k ) - G T k g k ^j 

+ a^Sofc + \\y k \\ 2 Rk - 2ylR k c k + \\c k \\ Rk - g£g k , 
0o = 0. (27) 

B. Reverse-Time Filter 

Here, we derive the Riccati equation, the filter state equa- 
tion, and the level shift scalar equation for the reverse-time 
filter over the time interval [k,t]. This comprises of solving 
for the infimum involving the function S2O) m dl4l . This is 
given by 

Bk(x) := inf S 2 (x k ,wi.)), (28) 

«)(.) ££ 2 [fc,*] 

where S^-) is defined in (fT2l . The corresponding reverse-time 
discrete-time HJB equation for the system (0 and the SQC 
(fT2l over the time interval [k, t] is given by 

V k (x,y : k ) = inf (14 +1 (a fc (a;) + D k w k ) 
iu ( .)e£ 2 [fc,t] 

|2 i ||„, ||2 i |2\ 



with 



V t (x) = 0. 



(30) 



As mentioned in IV-AI a first step towards obtaining an 
approximate solution to the HJB d29| i, we approximate the 
value function V k (-) as 



V k (x) w ||x - ijfellg + V'fc- 



(31) 



Applying the approximate solution (|3T| i to the HJB d29b we 

get 



\ x ~ ^IlL + V^fc = inf (||o!fe(^) + V'fe+i 

fc ti) ( .)6£ 2 [fc,i] fc + 1 

\z k \ 2 ). (32) 



+ IMI| fc + IMIk 



Solving for the infimum problem on the right hand side of 
(l32l using the method of completing the squares, we obtain 
the following optimum value for w k , denoted as w* 



Substituting w k 
terms we get 



from d33l l into d32l and rearranging 



--0fc = ||afc(a;)||n ^ 2x k D,a k (x) 

+ \\xk\\n + \\Vk -€k(x)\\%, - \k(x)\ 2 + V'fe+i, 

(34) 



where 



n = n fe+1 - u k+1 D k (Di u k+1 D k + Q k )*r>ln k+1 . (35) 



Here, (•)# denotes the Moore-Penrose pseudo-inverse. If the 
matrix (D k TL k D k + Q k ) is positive-definite, the pseudo- 
inverse in d35l l can be replaced by the normal matrix inverse. 
This is feasible and holds in most cases where a suitable 
solution exists for Yl k ; see e.g., JSJ. 

In order to obtain an approximate solution to the right hand 
side of (l34l l. we consider first order linearization of various 
nonlinear functions in (l34l l about the point x k as 



a k (x) w a k + A k (x - x k ), 
K k (x) w k k + G k (x - x k ), 
tik{x)~£k + C k {x-x k ), (36) 



where a k = a k (x k ),k k = n k (x k ), and £ k = Also, 
A k ,G k , and Cfc are the Jacobian matrices of the maps a k , K k , 
and respectively. 

Substituting d36l > into d34l l we obtain 



x 



2x 



I OA fc + Ai n + - G fc G fc i fe 



' iFfcl 
2^ 



- A T k m k + Cji?^ - | fc ) + G\k k 

\\ 2 AlnA k +n+2nA k +cl R k C k -GlG k 
[A T k m k + Qa k - ClR k {y k - &) - G^ fc ) 



10 



(Di U k+1 D k + Q fc ) n fc+1 (i fe - a k {x)). (33) 



"fc^fc + hk\\ Rk - 2y fc i?fcCfc + II&IIh* - KfeKk + Vfc+i- 

(37) 



Comparing terms on the left hand side and right hand side 
of (|37| |. we obtain the following recursion equations for the 



reverse-time filter corresponding to the time interval [k,t}\ 
Riccati equation over the time interval [k,t] : 
fife = A^flAk + C^RkCk — G^Gk, 
fit = 0. (38) 
Filter state equation over the time interval [k, t] : 
ik+i = Xk - flfc : A, 

x t = Xf (39) 
where 

A = (Aln(i k +i - &k) + ClR{y - £k) + GlGk). 
Level shift scalar equation over the time interval [k, t] : 

4>k+l = Ipk - \\Xk\\\T nAk+n _ 2nAh+C T RkCk _ G T ak 

+ 2xl(AlSla k + Qa k - C^R k (y k - &) - G T k k k ) 

- a k fla k - \\y k \\ 2 Rk + 2y k r R k £k - \\lkf Rk + K k k k , 
V>* = 0. (40) 

Hence, the robust smoothing algorithm will involve solv- 
ing for two sets of Riccati differential equations, filter state 
equations, and level shift scalar equations. These correspond 
to (f25t-(f27T> for the forward-time filter over the time interval 
[0, k] and (I3~8l-(l4"0b for the reverse-time filter over the time 
interval [k,t]. 

Consider the recursion equations (|25t-(|27t to have solutions 
x k ,H k ,<fik over the time interval [0, k], such that life ^ 
0, Qk ^ 0, Rk ^ 0. Also, if the recursion equations (l38l - 
d4Qb have solutions x k ,Tl k ,ijj k over the time interval [k,t], 
such that n^, ^ Q,Qk ^ 0, R k ^ 0, then the corresponding 
set-valued state smoother is given by 

X k [xo,y[o,t],d] = \x k e K" : ||x - Xfclln* + \\ x ~ %k\\n„ 

<d-4> k -ipk\. (41) 

Here, the level shift scalar equations for the forward-time 
filter in (|27| > and the reverse-time filter in (l40l are used to 
complete the definition of the set-valued state estimator (|4TT > 
for the robust smoothing algorithm and do not contribute to 
the Riccati difference equations or the filter state equations for 
either the forward-time or the reverse-time filters. 

VI. Conclusion and Future Works 

This paper derived a robust smoothing algorithm for a class 
of discrete-time uncertain nonlinear systems. It was assumed 
that the systems considered had uncertainties in modeling as 
well as exogenous noise sources. These uncertainties were 
described in terms of a sum quadratic constraint. The robust 
smoothing problem was solved by computing recursion equa- 
tions for a forward-time and a reverse-time filter. Both filters 
were formulated as set-valued state estimators which were 
further recast into corresponding optimal control problems. 
These optimal control problems were represented in terms 
of discrete-time Hamilton-Jacobi-Bellman equations. This for- 
mulation was straightforward in the case of the reverse-time 
filter. However, in the case of the forward-time filter, it was 



necessary to employ a forward-time formulation for the corre- 
sponding Hamilton-Jacobi-Bellman equation with reverse-time 
system dynamics, in order to ensure the forward-time recursion 
equations for the filter. The final solution to the smoothing 
problem consisted of computing approximate solutions to the 
underlying dynamic programming equations. This solution led 
to two sets of recursion equations and involved a Riccati 
difference equation, a filter state equation, and a level shift 
scalar equation for the forward-time as well as the reverse- 
time filters. 

As a note on future work, the robust smoothing algorithm 
presented in this paper can be applied to uncertain systems 
with unmodeled dynamics where it is needed to estimate 
outliers from noisy data. Such systems generally arise in the 
field of signal processing. 
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